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K strings in Yang-Mills theory can be considered as bound states of k elementary confining strings 
carrying one unit of colour flux. Current estimates of k-string tension Gic are very sensitive to 
the leading corrections due to quantum fluctuations of the string. In this study we address this 
problem by comparing Polyakov-Polyakov correlators in the fundamental representation {k= 1) 
with the corresponding ones with A: = 2 in the confining phase of a Z4 gauge theory in three 
dimensions. Highly efficient simulation techniques are available in this case. Although the k — I 
Polyakov-Polyakov correlator matches nicely with the expected bosonic string effects up to the 
Next-to-Leading-Order, the k ^2 Polyakov-Polyakov correlators show large deviations. This is 
an important source of potential systematic errors in the current estimates of a^. 
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1. Introduction 

Significant effort has been invested recently in studies of the flux tubes induced by color 
sources in higher representations of SU{N), built up of j copies of quarks in the fundamental 
representation. The long-distance properties of the flux tube should depend only on its N—ality 
k = j{modN), since all representations with the same k can be converted into each other by the 
emission of an appropriate number of soft gluons. As a consequence, the heavier strings of given 
N-ality k are expected to decay into the string with smallest string tension. The corresponding 
string is usually referred to as a ^-string. If its tension Ot for any allowed k satisfies the inequality 
(^ki+k2 < (^ki + <^k2^ the ^c— String is stable against decay into two strings of lower A/^-ality. 

Stable k strings are expected to belong to the antisymmetric representation with k quarks. This 
fact can be simply understood in terms of Casimir scaling, i.e. the hypothesis that the string tension 
for a given representation is proportional to the quadratic Casimir operator [|l|, ^ : within the set of 
aU representations of N—ality k the antisymmetric one corresponds to the minimum of the Casimir 
eigenvalues, suggesting 

^^,)_^ kiN-k) 
N-l 

Another competing hypothesis is the sine law: 

^(,.) _ ^ sin{k7i/N) 
" sin(7r/A^) ' 

which has been derived in the large A'^ limit of .yV = 2 supersymmetric SU (N) gauge theory softly 
broken to ^ = 1 [|3|], in the M theory description of ^ = 1 supersymmetric SU (N) gauge theory 
[0] and, more recently, in the AdS/CFT correspondence In some cases this formula is expected 
to be exact, while in others the calculated values of Ok turn out to be slightly smaller than ajf\ 

Lattice calculations in pure SU (N) gauge models for A/^ = 6 and A'^ = 4,5,6,8 [^|8|] in 
D = 3 + 1 point to the string tensions lying partway between the Casimir scaling and the sine 
law, however there is no complete consensus and some dedicated studies favour the sine formula 

©■ 

These calculations extract the string tensions Ok from Polyakov correlators at a finite temper- 
ature T assuming the free string prediction [p^ 



(^k = O \r / • (1-1) 



(^k = ^ -.-'J > (1-2) 



Ok{T) = Ok-^ + 0(r4), (1.3) 

D 

There is good evidence that this is a very accurate approximation \.ok= \ strings in all confining 
gauge theories irrespective of the gauge group once T is small enough [|6[ |ri|, 12], but there is no 



evidence that this can be also extended for > 1. Actually we show in this work that for ^ > 1 



Eq.(|1.3|) is not adequate and needs corrections. 

Our starting point was the observation that the effective string description, at least in the case 
/: = 1, is believed to be universal, that is to hold for all confining gauge theories. On the other 
hand in order to check the validity of Eq.( |1.3[ ) it is important to do a much more accurate finite 
volume study than any currently available for SU [N) ^—strings. Hence we decided to study a 
three-dimensional Z4 gauge theory, which is the simplest model supporting a ^ = 2 string. Using 
a duality transformation^ it is possible to map this model into the symmetric Ashkin-Teller (AT) 
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model, where very high precision can be achieved on large lattices, through a non-local cluster 
algorithm. The outcome of this analysis is uncontroversial: though at low temperature the tension 



of the fundamental string fits nicely with Eq.(|1.3[) even at Next-to-Leading-Order (NLO), we find 



a clear mismatch for the 2-string in the same temperature range, suggesting the need of corrections 



toEq.(1.3) 



Theoretically, there is something one can say about the origin of these corrections. strings 
can be viewed as bound states of k fundamental strings. Accordingly, we expect that besides the 



mechanical vibration modes whose quantum contributions yield exactly Eq.(1.3), there should be 



breathing modes related to the internal degrees of freedom of the k constituent strings. 

2. The effective string model 

The infrared description of any confining gauge theory is well described by an effective string 
model. A particularly simple string action is the Nambu-Goto one, where the correlation function 
of two Polyakov loops at a temperature T = l/L and at a distance R can be calculated at the NLO 



[13, 14], yielding 



{P{0)pHR)) ^^^2 ■,r^-;E = 2E,-El (2.1) 



where T] is the Dedekind eta function and E„{x) are the Eisenstein functions (see e.g. [14] for 
detailed definitions). There is strong evidence that at this order this formula is universal, i.e. it 
holds in whatever confining gauge theory [15, T^, 17]. 



3. Algorithm 

We work directly in the dual form of the 3dZi4 gauge model, i.e. a symmetric Ashkin-Teller 
(AT) model. It is described in terms of two coupled, ferromagnetic, Ising systems through the 
two-parameter action 

5 = - ^ j8 {OjcGy + Z^Ty) + a O^OyX^Ty, (3. 1) 

to') 

where Ox and are the Ising variables (a^, = ±1). The global Z4 symmetry of the action is 
generated by the transformation a — > — T, T ^ a. An independent Z2 symmetry is generated by 
the transformation a T, which is related to the charge conjugation of the corresponding dual 
model. 

The great advantage of studying the 2>d AT model instead of the original 3<i Z4 gauge theory is 



that a non-local cluster updating algorithm [18] can be used. Of course, in the AT model we have 
two spin variables, so it is necessary to extend the original method. The idea is the following: we 
take the two site variables a and T as if they belonged to two distinct lattices Ra and R^; we freeze 
the lattice 7?^ and apply the Swendsen-Wang algorithm to the variables T so obtaining a new lattice 
R'^. At this point we freeze the lattice R'^ and we update the variables a, and so on. 

It is possible to show that in this statistical system the expectation value of the Wilson loop of 
the gauge model (Wj) gauge is given by: 

7* 



(yVy) gauge — , 
^AT 
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where Zat is the partition function of the AT model and Z^j is that modified by a suitable twist of 
the couplings. In particular, in order to determine the (Wy) gauge in the fundamental representation 
= 1) it suffices to flip the couplings of a or T just in the links orthogonal to an arbitrary surface £ 
encircled by 7 = 5 £. Similarly, flipping the signs of both couplings of a and T in the same surface, 
we get the Wilson loop in the k = 2 representation. 

Actually, since we have a model written in terms of Ising variables we can use a very powerful 
method to estimate the expectation value of the Wilson loop based only on the topological linking 
properties of the Fortuin-Kasteleyn (FK) clusters [|T9|]. For each FK configuration one looks for 
paths in the clusters and then applies the following rule: 

• = 1 if there is no path hnked with the loop 7 or if the winding number modulo 2 is zero, 

• Wy = otherwise. 

Note that using such a definition, the value of Wy does not change if dangling ends or bridges 
between closed path are added or removed [|20[]. 

In the AT model we apply the above rules separately for the two variables a and T, therefore 
we have, for each configuration, a value and a value W^, each of them corresponds to the value 
of the Wilson loop in the fundamental representation; the product W^W^ corresponds to the value 
of Wilson loop in the k = l representation. 

The same ideas do apply to measure the Polyakov-Polyakov correlator (P(O)P^ (/?)). 



4. Numerical results 



In this paper we discuss data obtained by simulating the AT model in one point of the parameter 
space (a = 0.05, j3 = 0.207), where our best estimate for string tensions are Oi = 0.02084(5) and 
(72 = 0.0323(5). In this case, the tension ratio is G2/01 ~ 1.55, while the values predicted by the 
sine or Casimir scaling for N = 4 and ^ = 2 are \/2 ~ 1.41 and 4/3, respectively. 

We performed 10^ measurements for both Oi and 02 on a lattice 64^ x L. The value of L has 
been determined in such a way that T ^r^Tc/l using the relation Tc/\/g ^ 1.1. 

In order to check whether the functional relation (2J_) is an accurate description of the Polyakov- 
Polyakov correlator two conditions are needed: the fitting parameters (in particular a^) should be 
stable at large R and the estimated value of Ot should not depend on T since Eq.(^?T|) should ac- 
count for the T dependence. 

In Fig. [l]a it is possible to see the Gf value^ fitted by the Leading-Order (LO) approximation 
(i.e. neglecting the E term in (2J_)) in a range [/?„„•„, 18] (where is the value that appears on the 
x-axes); albeit for large values of it seems to appear a plateau, there are three different values 
of Of, therefore the functional relation (FR) used is not sufficiently accurate to determine it. 

Fig. |I|b is similar to Fig. [l|a, but in this case we have used the same FR to determine the value 
of O//; also in this case we can apply the previous considerations and the FR is not correct. 

Then, in Fig. we fitted the value of Of to two different FRs, one obtained by a LO approx- 
imation and the other one by NLO; the value of Of increases significantly and the plateau appears 
for a smaller value of 



In this Section O/ = tJi and Off = 02- 
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Figure 1: oy and O// vs Rmin for three different T at LO 



In order to shed some light on the difficulty in determining a string tension that is free of 
systematic errors, in Fig. || we plot the values of x^/d.o.f. for fits which appear in Fig. ^ note that 
for Rmin > 1 1 there is no difference in x^/d.o.f. between LO and NLO, even if the values of Of in 
the same range (see Fig. ^ are different. It is interesting to note that in NLO case also the value of 
X^/d.o.f. shows a plateau for Rmin > 9, therefore this FR better describes what happens at small 
values of 7?. 

Fig. ^ is similar to Fig. [l^ but now we use FRs obtained with a NLO approximation; now the 
FR is actually correct: the value of Of is stable also for small values of Rmin and there is only one 
value of it for three different values of L (i.e. different values of T). 

Finally, in Fig. ^ and Fig. ^ we interpolate the value of Off with NLO approximation with 
two different hypotheses: in Fig. |5]a the two strings are stuck together and fluctuate as a single 
string (we put D=3 in (p!l|)); in Fig. ^ the two strings are assumed to fluctuate independently 
(this is equivalent to putting D=4 in (|2.1D). It is clear that in this case the FR does not describe 
accurately Of j: in both cases a plateau does not appear for a// and there are different values for 
different values of T. This indicates that the current estimates of Ok are affected by systematic 
errors. 
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Figure 2: Gf vs for LO and NLO 
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Figure 3: ;t^^ /d.o.f vs for LO and NLO 
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Figure 4: (7/ vs for three different T at NLO 



5. Conclusions 

The functional relations which are very accurate to fit data related to fundamental string tension 
are not adequate to describe 2-string tension. 

An accurate estimate of the ^-string tension is rather problematic: from a numerical point of 
view our analysis shows that a blind application of the usual formulas for the Polyakov-Polyakov 
correlators introduces strong systematic errors. From a theoretical point of view it is necessary to 
face the problem of determining the correct functional relation of Polyakov-Polyakov correlator for 
higher representations, which is presently unknown. 
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